UNCLASSIFIED 


AD  NUMBER 

AD845078 

NEW  LIMITATION  CHANGE 
TO 

Approved  for  public  release^  distribution 
unlimited 


EROM 

Distribution  authorized  to  U.S.  Gov't, 
agencies  and  their  contractors; 
Administrative/Operational  Use;  DEC  1968. 
Other  requests  shall  be  referred  to  Army 
Missile  Command^  Redstone  Arsenal^  AL. 


AUTHORITY 

USAMC  Itr,  23  Aug  1971 


THIS  PAGE  IS  UNCLASSIEIED 


Teclmic<ti  Beport  8-lBi 


DB^SLOPBSKT  OF  A  BALL  2®SNTATK)N  l^St 
FCm  SASIPI^S  Finns  THICKHB^ 

by 

WEUam  M.  Sigmon 
Ciiarles  H.  Parr. 

Albert  J.  Jgnatowski 


December  1868 


■i- 


Ue  8,  ARMY  MISSJLE  COMMAND 
HEBSTOKE  ABSEKAL  ALABAMA  35803 


J 


Contracts 

DAAH0i-67-C-0655 

DAAH01-66-C-08S2 


REDSTONE  RESEARCH  L^ORATORIES 

HUNTSVILLE,  ALABAMA  35897 


^0^  ^  fe^  Ttmd^  In  acc^^acee 

V%i%  irt^emasts  rs- 

$;^v^  ^^42^  goTSmsio^ 

-•^^t  ^tac  i3j9  Da¬ 
ises#  Eagtllng  5, 


‘I  i2  to  es^Knrt  ss^  eacSi 

tj^SfSdi®  to  Cc£sl^-^|gi2ia$3ts  or  ^^Iga  mtiffials  may  bo 
133^3  o?%  'STiib'p^i^  a^oyni  c!^ 


liK^'iiSSr 


gB2^7s  4ERK1  Of 
rjf  S3I?55SS| 
ESRin  £j| 

ft6.,il3l  _ _ I 


Ife^apliasat  o!  Army 

tJ,  S,  Army  MissUe  Ccmimajig 
^sfcao  Aresi^aJ,  Alaiama  3^0 


DISCLAIMER 


to  fiaSiaags  in  report  are  not  to  be  constrnod  as  an  o!a~ 
cial  D^sariment  of  the  Army  position  naless  so  designated  by 
other  antJaorisefi  docna^ats. 


il 

tfc 


DISPO^^TB^N  MSTRUCTIDNS 

■t.coy  ih\Q.  :re^;-j4'i  ,,htT'  no  ior»gQr  rteeocC',.  Do  not  return  it  tc 
CjOi^'Aolor-. 


December  1968 


Tecimical  Report  S-182 


DEVELOPMENT  OF  A  BALL  ]NDENTATK>N  TEST 
FOR  SAMPLES  OF  FINirE  THICKNESS 

! 

i 


by 


William  M.  Sigmon 
Charles  H.  Parr 
Albert  J.  Ignatowski 


STAIEi£Z::T  #2  UK-CLASSIFIED 


This  docunent  Is  subject  to  special  export  controls  and  eadi 
transmittal  to  foreign  governments  or  foreign  nationals  Bay  be 

aade  only  fith  priar  eypxvftl  ■■■■■■■■■■■-■■  ■■  ■ 


U.  S.  ARMY  MISSILE  COMMAND 
Redstone  Arsenal,  Alabama  35809 


Contracts 

DAAH01-67-C-0655 

DAAH01-68-C-0632 


Distribatioo  Ltnited 


ROHM  AND  HAAS  COMPANY 

REDSTONE  RESEARCH  LABORATORIES 
HUNTSVILLE,  ALABAMA  35807 


FOREWORD 


The  work  described  in  this  report  was  performed  onder 
Contracts  DAAH01-e?-C-Ofc55  and  DAAH01-68‘  C-0632  for  exploratory 
development  of  propellants  for  missiles  and  rockets  under  the  techni¬ 
cal  cognisance  of  the  Solid  Propellant  Chemistry  Branch.  Army 
Propulsion  Laboratory  and  Center,  Research  and  Development 
Directorate.  O.  S.  Army  Missile  Command. 

The  propellant  chemist  olton  rc^cests  the  mechanical  testing 
laboratory  to  supply  some  form  of  quantitative  laboratory  meastire 
of  the  mechanical  behavior  of  small  quantities  of  low  modulus  propel¬ 
lant  binder  gumstocks.  Frequcrtly  these  materials  are  ill-suited 
for  the  usual  gamut  of  mechanical  property  tests  developed  for 
finished  propellants.  Mechanical  property  data  supplied  at  an  early 
phase  of  propellant  binder  development  can  be  useful  in  determining 
proper  curing  systems  and  conditions,  effects  of  plasticisers  or  other 
modifiers,  effects  of  variations  of  binder  synthesis  or  processing, 
and  can  give  indications  of  glass  transition  temperature,  strain  rate 
dependency,  and  other  factors  which  are  important  in  the  finished 
propellant.  The  test  developed  here  is  useful  in  these  applications 
and  deviates  from  previous  hardness,  indentation,  or  flexure  tests 
in  that  the  accompanying  data  analysis  includes  the  effect  of  finite 
sample  dimensions. 

The  work  described  here  grew  out  of  work  first  initiated 
under  Contract  DA-01-021  AMC- 1 1536(Z). 


ABSTRACT 


A  spherical  indentatioc  test  for  sxztall  quantities  of  low  modtdus 
materials  is  described  and  experimentally  validated.  The  data  analy¬ 
sis  takes  into  accoon*^  finite  sample  fliickness.  Comparison  with  the 
Hertz  contact  theory  is  made,  and  extension  to  viscoelastic  charac¬ 
terization  is  discussed. 
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Section  I.  INTRODUCTION 

A  ball  indentation  test  has  been  developed  to  meet  the  need  for 
a  simple,  practical  method  of  measuring  the  shear  modulus  of  labora¬ 
tory  q»’a'C!titie;-  of  rubbery  materials.  The  indentation  test  offers  a 
number  of  possible  advantages:  sample  sice  and  shape  are  noncritical 
and  only  a  rvmiimum  of  preparation  is  needed;  the  quantities  of  sample 
material  required  are  relatively  small;  die  test  device  is  simple  and 
easy  to  operate-  the  method  is  non-destructive  in  character,  thus 
permitting  repeat^'d  tests  on  the  same  sample. 

Many  of  the  previous  applications  of  the  indentation  test  were 
limited  by  lack  of  a  well-defined  relationship  among  applied  load, 
indentor  geometry,  indentation  depth,  and  the  material  shear  modulus, 
particularly  for  thin  test  samples.  Most  work  in  the  past  made  use 
of  either  the  classical  Hert*.  contact  theory  or  of  some  modificatian 
of  it  based  upon  empirical  te.^t  data  [1-4].  The  Hertz  solution,  dis¬ 
cussed  in  many  texts  on  elash  uty  (e.  g. ,  Timoshenko  and  Goodier 
[5]),  assumes  the  indented  solid  to  be  a  semi- infinite  layer,  i.e. , 
of  infinite  depth  below  the  defoii.-ced  surface.  Therefore,  in  the 
design  of  a  practical  teat,  questions  arise  as  to  (1)  the  limits  of  finite 
specimen  size  and  indentation  for  which  the  Hertz  theory  will  apply 
with  sufficient  accuracy,  and  (2)  other  relationships  which  may  be 
used  for  those  cases  where  the  Hertz  theory  will  not  provide  results 
of  the  required  accuracy. 

Hitherto,  mainly  because  no  precue  information  was  available, 
the  classical  theory  was  as,?umed  to  hold  f  yr  an  indentation  not  ex¬ 
ceeding  10  percent  of  the  specimen  thicknesj'-  [1,2,4].  Gent  [2],  for 
example,  has  shown  that  the  Hertz  theory  is  qjite  accurate  for  de¬ 
termining  Young's  modules  of  vulcanized  rubbers  when  using  small 
indentations,  up  to  one-tenth  the  sample  thickncf^s.  Waters  [4]  con¬ 
ducted  a  series  of  tests,  using  spherical  indei  tor?  on  sheets  of  vul¬ 
canized  rubber  of  various  thicknesses,  and  fouKo  tiv^t  moduli  com¬ 
puted  by  the  Hertz  theory  were  in  error  by  less  than  5  percent  when 
the  sheet  thickness  was  more  than  8  times  t-he  indsntor  contact 
area  radius.  With  thinner  sheets,  much  larger  errors  were  experi¬ 
enced,  exceeding  30  percent  for  thicknesses  less  than  times  the 
contact  radius.  Waters'  approach  to  the  problem  war  to  ..oodify  the 
Hertz  expression  by  an  experimentally  derived  function  involving  the 
ratio  of  thickness  to  contact  area  and  including  a  consta;..t  it  Allow 
for  the  surface  conditions  (lubricated  and  dry). 

None  of  these  methods  are  particularly  applicable  tc  the 
problem  of  measuring  the  shear  modulus  of  "soft"  polymers  'thoi  e 


with  moduli  less  than  lO’  psi)  in  the  form  of  relatively  thin  sheets. 

In  order  to  have  indentation  depths  and  forces  sufficiently  large  to 
be  recorded  accurately  by  conventional  instruments,  the  specimen 
must  be  indented  to  a  considerable  fraction  of  its  total  thickness, 
tixerefore  deviating  ir*arkedly  from  the  Hertz  model.  Empirical  for¬ 
mulas  are  of  doubtful  value,  considering  the  wide  range  ol  materials 
encounteied.  However,  Lebedev  anc  Ufliand  [6]  have  dealt  with  the 
problem  of  pressing  a  rigid  axisymmetric  punch  of  any  profile  into 
an  elastic  layer  of  finite  thickness.  Their  solution  offers  the  possi¬ 
bility  of  obtaining  a  firm  theoretical  relationship  among  load,  in¬ 
dentation,  and  shear  modulus  for  tests  cf  rubbery  materials  in  readily 
available  forms.  Parr  [7,  8]  has  applied  the  Lebendev- Ufliand  solution 
to  the  specialized  case  of  a  spherical  indentor,  with  particular  em¬ 
phasis  on  obtaining  a  relationship  in  a  form  most  useful  to  the  test 
laboratory.  The  results  of  this  analysis  led  to  the  development  of 
the  ball  indentation  test  method  described  herein. 
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Section  II.  ANALYSIS  OF  THE  INDENTATION  TEST 
I .  Elastic  Materials 

The  ball  indentation  test  may  be  considered  as  a  special 
case  of  the  general  problem  of  a  punch  penetrating  a  solid.  The 
geometry  is  shown  in  Figure  1.  The  object  of  the  analysis  is  to  find 
a  solution  relating  the  unknown  shear  modulus  of  the  test  sample  to 
the  known  or  measurable  quantities  of  indentor  radius  R,  applied  load 
P,  and  penetration  Wq,  for  any  thickness,  h,  of  sample  material.  The 
materials  to  be  tested  by  this  method  will  have  elastic  moduli  several 
orders  of  magnitude  smaller  than  those  of  steel  or  aluminum,  of 
which  the  test  device  is  constructed;  therefore,  for  practical  purposes 
the  indentor  and  base  may  be  considered  infinitely  rigid. 

The  solution  to  the  problem  is  outlined  in  the  Appendix,  Briefly 
stated,  the  solution  consists  of  using  the  method  developed  by  Licbedev 
and  UfLiand  to  express  displacements  and  stressL.s  in  terms  of  one 
auxiliary  function,  which  represents  the  solution  of  a  Fredholm 
integral  equation.  The  general  solution  is  reduced  to  the  specific 
case  of  a  spherical  indentor,  and  equations  are  obtained  for  applied 
load  and  indentor  penetration  in  terms  of  the  auxiliary  function.  The 
expressions  thus  obtained  are  too  complex  to  be  evaluated  exactly, 

»  therefore,  numerical  methods  are  employed  to  compute  values  of 

the  auxiliary  function.  These  values,  in  turn,  are  used  to  solve  for 
load  and  penetration  by  quadrature. 

The  results  are  presented  in  Figures  2  and  3.  These  sets  of 
curves  are  taken  from  the  paper  of  Ignatowski  and  Parr  [8]  in  which 
the  analysis  originally  appeared.  In  Figure  2  the  dimensionless  load 
(l-v)P/uR^  is  plotted  versus  dimensionless  indentation  Wq/R  for 
various  values  of  the  dimensionless  thickness  ratio  h/R  (v  is  Poisson's 
ratin  and  fj.  the  shear  modulus).  Figure  3  presents  the  same  infor¬ 
mal  -n  in  a  different  manner,  with  indentation  as  the  abscissa  and 
curves  for  various  values  of  thickness  ratio.  Note  that  the  elastic 
constants  *  '.he  material  appear  only  in  the  dimensionless  load;  all 
computatic,-8  are  independent  of  the  elastic  constants.  This  situation 
results  from  the  boundary  conditions  imposed,  in  which  frictionless 
contact  between  the  test  sample  and  base  or  indentor  is  assumed. 

Dimensionless  quantities  were  chosen  for  the  calculations 
and  Figures  2  and  3  for  the  sake  of  convenience  and  generality,  but 
»  these  quantities  may  not  neces.iarily  be  the  most  convenient  for 

actual  test  data  reduction.  The  curves  can  easily  be  modified  to 
apply  to  any  particular  laboratory  setup.  With  the  indentor  radius 
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FIGURE  1.  GEOMETRY  OF  INDENTATION  PROBLEM 
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known,  and  Poisson's  ratio  known  or  assumed  (  v  =  0.  5  is  a  very 
good  approximation  for  rubbery  materials),  curves  of  P/ji  versus 
Wq  for  various  values  of  h  (from  Figure  2),  or  versus  h  for  various 
values  of  Wq  (from  Figure  3),  are  obtained.  The  shear  modulus 
can  then  be  determined  from  laboratory  measurements  of  P,  h, 
and  Wq.  a  difierent  set  of  curves  is  required  for  each  indentor 
radius  and  Poisson's  ratio,  but  these  are  available  simply  by  changing 
the  scales  on  the  basic  curves  in  Figures  2  and  3. 

The  curve  of  the  classical  Hertz  contact  problem  solution 
is  included  in  Figure  2  for  comparison  with  the  Lebedev- Ufliand 
solution.  Note  that  the  two  essentially  agree  for  specimen  thick¬ 
nesses  in  excess  of  ten  indentor  radii.  However,  for  thinner  speci¬ 
mens,  drastic  differences  between  the  solutions  are  seen.  The  Hertz 
theory,  developed  for  a  semi-infinite  layer,  does  not  consider  the 
effects  of  finite  thickness  which  become  prominent  for  relatively  thin 
layers.  For  example,  for  a  specimen  thickness  of  one  indentor 
radius  (h/R  =  1)  the  Hertz  theory  gives  a  shear  modulus  value  some 
35  percent  higher  than  that  given  by  the  Lebedev- Ufliand  curves. 
Additional  comparisons  will  further  define  the  practical  limits  of  the 
classical  theory  and  illustrate  the  hazards  of  its  misuse  with  thin 
sections . 

2.  Viscoelastic  Materials 

Strictly  speaking,  the  foregoing  discussion  and  the  analysis 
given  in  the  Appendix  pertain  only  to  elastic  materials.  However, 
since  many,  perhaps  most,  materials  of  interest  exhibit  viscoelastic 
behavior,  it  becomes  important  to  know  how  the  indentation  test  may 
be  applied  in  measurement  of  viscoelastic  properties. 

Through  the  processes  of  stress  relaxation,  certain  visco¬ 
elastic  materials,  subjected  to  a  constant  strain,  will  eventually 
reach  a  condition  of  equilibrum  for  which  the  stress  remains  constant 
with  time.  Moduli  corresponding  to  this  condition  are  known  as 
equilibrium,  long-time,  or  "rubbery"  moduli,  and  may  be  computed 
from  indentation  test  data  and  the  results  of  the  elastic  analysis, 
just  as  with  an  elastic  material.  In  other  words,  the  material  i.*? 
assumed  to  have  reached  a  state  where  the  shear  modulus  is  indepen¬ 
dent  of  time  as  it  is  in  an  elastic  material.  The  only  requirement 
is  to  wait  a  sufficient  length  of  time  before  taking  readings  to  in¬ 
sure  that  equilibrium  has  been  reached. 

For  shorter  ranges  of  time  when  equilibrium  stress-strain 
conditions  have  not  be  achieved,  or  for  tests  where  the  purpose  is 


measurement  of  short-time  effects  such  as  relaxation  moduli,  the 
indentation  test  and  its  elastic  solution  also  appear  useful.  However, 
the  problem  of  application  becomes  considerably  more  difficult. 

The  information  in  Figure  2  maybe  considered  as  described 
by 


P 


H'^o) 


(1) 


where  f(w(j)  is  a  function  known  only  graphically.  Lee  and  Radok 
[9]  have  shown  that  the  viscoelastic  counterpart  of  tne  Hertz  con¬ 
tact  theory  cam  be  deduced  from  the  elastic  olution.  By  analogy, 
we  may  say  that 

P(t)  =  f  H{t-T)  dT  (2) 

where  |i.(t)  is  the  shear  relaxation  modulus,  and  f(wo)  is  time  depen¬ 
dent  because  Wq  is  time  dependent.  Equation  (2)  may  be  solved  for 
u(t)  if  the  time  dependency  of  f(wo)  is  known.  This,  however,  is  not 
easily  done  even  if  Wo(t)  is  known  because  f(wo)  is  a  complex  function 
of  Wo(t), 

One  useful  case,  for  which  there  is  a  relatively  simple  solution, 
is  that  of  a  step  function  input.  If  Wo(t)  is  applied  as  a  step  function, 
or  instantaneous,  penetration,  then  f(wo)  is  also  a  step  function  and 
(2)  reduces  to 

P(t)  =  p(t)f(wo)  .  (3) 

The  dimensionless  load  in  Figure  2  can  then  be  interpreted  as 

(l-v)P(t) 

RV(t) 

Knowledge  of  the  load-time  history,  P(t),  permits  calculation  of  the 
shear  relaxation  modulus  (i.(t). 

In  actual  practice  a  step  function  is,  of  course,  not  achieved. 
Some  finite  time,  tj,  is  required  for  the  indentor  to  reach  the  de¬ 
sired  penetration  depth.  Nevertheless,  for  many  materials  it  may 
be  possible,  with  suitable  high-speed  test  devices,  to  reduce  tj  to 


8 


values  much  less  than  the  times  of  interest  and  thus  be  able  to  consider 
Wo(t)  a  step  function  for  all  practical  purposes.  Following  the  saune 
approximation  often  used  in  tensile  stress  relaxation  testing  [1C],  it 
appears  reasonable  to  assume  a  step  function  input  for  times  greater 
than  10  t|. 
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Section  III.  DESCRIPTION  OF  EXPERLK-IENTAL 
EQUIPMENT  AND  PROCEDURES 


1 .  Test  Device 

A  simple  prototype  test  device  was  initially  built  to  obtain 
a:,  approximate  comparison  between  the  theoretical  analysis  and 
experimental  data.  The  results  were  encouraging,  and  so  a  more 
refined  design  model  was  constructed  according  to  the  sketch  of 
Figure  4.  The  ball  indentation  test  device  shown  in  Figure  4  some¬ 
what  resembles  two  interleaved  C-clamps.  The  upper  "C"  supports 
the  test  specimen  and  is  designed  tc  be  suspended  from  a  load  cell 
by  means  of  a  pull  rod  with  universal  joints .  It  is  counterweighted 
and  balanced  so  as  to  hang  freely  with  the  pull  rod  centerline 
vertical.  The  spherical  indentor  is  mounted  in  the  lower  "C," 
which  is  attached  by  another  pull  rod  to  the  movable  crosshead  of  a 
universal  tensile  testing  machine.  With  the  test  specimen  and  in¬ 
dentor  located  between  the  overlapping  arms  of  the  "C  ”s ,  the  in¬ 
dentor  is  pressed  into  the  specimen  when  the  lower  ”C"  is  given  a 
downward  displacement.  Scribe  marks  along  the  centerline  of  both 
parts  facilitate  alignment  during  assembly  and  operation. 

The  test  device  was  originally  designed  to  be  used  with  Instron 
universal  test  machines,  but  is  easily  adaptable  to  many  other  loading 
arrangements.  For  example,  a  fixed  support  and  a  dead  weight  pro¬ 
duce  a  creep  tester.  The  small  size  of  the  device  (approximately 
5X3X2  inches  overall,  not  including  pull  rods)  enables  it  to  be 
used  in  most  available  temperature  and  environmental  chambers. 

The  indentors  are  commercially  available,  ground  emd  polished, 
tungsten  carbide  balls,  cemented  to  insert  pins.  The  indentors  are 
readily  interchanged  and  are  presently  used  in  a  range  of  sizes. 

Table  I  gives  the  sizes  which  are  available.  Maocimum  test  specimen 
dimensions  are  1. 5  inches  square  by  0.  75  inch  thick.  The  test  speci¬ 
men  base  has  been  coated  v/ith  Teflon®-  to  reduce  friction. 

2.  Test  Procedures 


The  only  test  specimen  preparation  usually  required  is  to  cut 
the  specimen  to  a  size  which  will  fit  the  test  device.  Very  rough 
cuts  are  acceptable  for  the  edges,  with  no  restrictions  on  squareness. 
The  contact  surfaces,  between  the  specimen  and  the  base  and  indentor. 


^Trademark  of  E.  I.  DuPont  and  Company,  Wilmington,  Del, 
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Table  I.  Indentation  Test  Data  for  Hysol  8705  Specimen 


For  all  test  conditions; 

Specimen  dimension  approx.  1, 5  in.  sq. 

Specimen  thickness  =  h  =  0.  504  in. 

Poisson's  ratio  =  v  =0.46  from  Mfg's  specs 
Penetration  rate  =  2,0  in. /min 
Temperature  77  T 

Load  readings  taken  after  2  min  relaxation  at  constant 
indentation 

Specimen  surfaces  dusted  with  talcum  powder 


Indentor 
Radius  > 

R 

(in. ) 

h/R 

Indentation 
Depth,  Wq 
(in.) 

Wq/R 

iw 

Load,  P 
(Note  2) 
(lb) 

Computed 
Shear 
Modulus,  p 
(psi) 

0.0625 

3.05 

0.025 

mm 

0.69 

1.05 

212 

0.038 

IM 

1.25 

1.92 

212 

0.050 

0.80 

1.92 

2.97 

214 

0.0938 

5.37 

0.038 

0.40 

0.  71 

2.  7.3 

235 

0.056 

0.60 

1.30 

4.  74 

222 

0.075 

0.80 

1.97 

7.45 

231 

0,094 

i.OO 

2,  72 

9.92 

223 

0.  1250 

4.02 

0.050 

0,40 

0.  74 

4.78 

224 

0.075 

0.  60 

1. 34 

8.80 

227 

0.  1562 

3.22 

0.031 

0.20 

0.26 

2.45 

208 

0.062 

0.40 

0.  76 

7.  19 

209 

0.  1875 

2.69 

0.038 

0.20 

_ 

0.  26 

3,  56 

206 

0.2188 

2.30 

0.044 

0.20 

0.27 

5.  70 

234 

0.  2500 

2.01 

0,050 

0.20 

0.28 

7.  28 

224 

Mean  Value 

= 

220 

Standard  Deviation  = 

10 

Note  1  -  Dimensionless  load  values  from  curves, 
Note  2  -  AveraRe  reading  from  2  runs. 

Figure  3. 
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should  be  reasonably  flat,  smooth,  and  parallel.  Allowable  deviations 
depend  upon  the  particular  test  circumstances  regarding  desired  accu¬ 
racy,  ease  of  operations,  and  specimen  properties.  Care  must  be 
taken  about  non-parallelism,  however,  because  non-parallel  surfaces 
will  introduce  horizontal  force  components  which  may  cause  appreciable 
errors,  and,  in  extreme  cases,  even  cause  the  indentor  to  slide 
across  the  specimen  surface. 

It  is  also  suggested  that  the  specimen  surfaces  be  dusted  with 
talcum  powder  or  lubricated  in  some  other  fashion.  The  effects  of 
surface  conditions  have  not  yet  been  investigated  fully,  but  pre¬ 
liminary  checks  indicate  more  consistent  results  with  lubricated  sur¬ 
faces.  However,  the  possibility  of  contamination  of  the  test  sample 
must  be  considered  in  the  choice  of  a  lubricant,  and  if  any  doubts 
exist  lubrication  should  be  omitted  entirely.  The  effect  upon  the  test 
results  would  probably  be  minor. 

Other  limitations  exist  regarding  the  minimum  test  specimen 
dimensions.  In  the  theoretical  analysis  the  specimen  was  assumed 
to  be  infinite  in  extent,  though  of  finite  thickness.  No  analysis  has 
yet  been  performed  for  a  specimen  of  finite  dimensions,  but  test  re¬ 
sults  have  yielded  some  useful  guidelines.  Apparently  specimen  di¬ 
mensions  are  not  at  all  critical  so  long  as  the  length  and  width  are 
at  least  1. 5  times  the  thickness  and  6  times  the  radius  of  the  largest 
indentor  to  be  used.  Reasons  lor  choosing  these  values  will  be  dis¬ 
cussed  in  a  later  section. 

The  test  procedure  consists  of  displacing  the  lower  "C" 
fixture  and  recording  the  applied  load  and  indentation  depth  as  functions 
of  time.  Indentation  depths  may  be  determined  from  the  test  machine 
indicator  dials,  a  separate  displacement  transducer,  or  by  reading 
elapsed  time  from  the  load-time  record  and  multiplying  by  cross - 
head  speed  (if  known  and  constant). 

The  fractional  or  percentage  error  in  determining  indentation 
depth  may  be  minimized  by  going  to  as  great  a  depth  as  possible, 
subject  to  the  limitations  imposed  by  specimen  thickness  and  allowable 
loads.  Penetration  must  not  exceed  the  spherical  radius,  because 
the  equation  used  for  the  contact  surface  contour  in  the  analysis  would 
then  no  longer  apply.  Smaller  indentor  sizes  result  in  lower  loads 
at  a  given  indentation,  but  are  limited  in  penetration  depth;  therefore, 
some  compromise  may  be  necessary. 


Section  IV.  EXPERIMENTAL  RESULTS  AND  DISCUSSION 


I .  Evaluation  of  Test  Device  and  Method 

A  number  of  different  materials  have  been  tested  with  the 
bail  indentation  test  ’evice,  but  the  results  of  only  a  few  will  be 
presented  herein  for  the  purpose  of  illustration.  Two  test  series 
of  interest  are  those  of  Hysol®^  and  Solithane®^  113  polyurethane 
rubbers.  These  materials  were  used  to  evaluate*  the  test  device  and 
method  because  their  properties  were  well  know'n  from  manufacturer's 
data  or  other  test  results.'*  Comparison  of  the  shear  moduli  obtained 
from  the  indentation  test  with  those  obtained  by  other  means  provided 
a  check  on  the  validity  of  the  analysis  and  the  curves  in  Figures  2  and 
3. 

Specimens  of  both  materials  were  tested  in  the  same  manner. 
The  indentor  was  pressed  at  a  constant  rate  into  the  specimen  to  a 
predetermined  depth,  and  held  constant  at  this  depth  for  two  minutes. 
The  curves  of  Figure  3  were  entered  with  the  values  of  Wq,  R,  and 
h  used,  along  with  the  load,  P,  read  after  two  minutes,  and  values 
of  the  shear  modulus  were  obtained.  This  modulus  was  considered 
the  equilibrium  shear  modulus,  since  an  essentially  constant  value 
of  load  was  reached  v/ell  within  the  two  minute  waiting  period. 

Data  from  tests  on  the  Hysol  8705  specimen  are  summarized 
in  Table  I.  Note  that  the  same  specimen  was  tested  with  a  number  of 
different  size  indentors  and  various  indentation  depths.  No  significant 
differences  in  results  were  found  for  any  of  the  different  indentors 
or  depths.  This  was  an  expected  result,  for  there  is  no  indication  in 
the  analysis  that  the  test  setup  parameters  should  affect  the  computed 
shear  moduli.  The  mean  value  of  shear  modulus  and  its  standard 
deviation,  calculated  from  the  test  data,  are  given  in  Table  I  as  220 
psi  and  10  psi,  respectively.  A  shear  modulus  was  also  computed 
from  Young's  modulus  and  Poisson's  ratio  according  to  the  familiar 
equation 

^Trademark  of  Hysol  Corporation,  Clean,  N.  Y. 

^Trademark  of  Thiokol  Chemical  Corp.  ,  Trenton,  N.  J. 

*The  specific  Solithane  sample  tested  was  provided  through  the 
courtesy  of  Professor  W,  G.  Knauss,  and  was  one  of  the  materials 
thoroughly  characterized  as  part  of  the  development  of  a  standard 
crosslinked  polymer  program  conducted  at  the  California  Institute  of 
Technology  under  U.  S.  Air  Force  sponsorship. 
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where  Young's  modulus,  E,  was  obtained  by  conventional  tensile  tests 
on  specimens  cut  from  the  same  sheet.  The  tensile  test  data  ^rve  an 
average  E  =  550  psi,  which  in  turn  gives  p  =  190  psi.  Agreement 
between  the  two  values  of  shear  modulus  is  considered  reasoncoiy  good. 


Data  from  tests  of  the  Solithane  1 13  specimen  are  presented 
in  Table  II.  Again,  a  number  of  various  size  indentors  and  inden¬ 
tation  depths  were  used,  along  with  different  penetration  rates.  As 
before,  no  significant  differences  were  noted  in  the  values  of  shear 
modulus  found  for  any  of  the  different  test  conditions.  The  shear 
modulus  mean  value  and  standard  deviation  were  calculated  to  be 
197  psi  and  16  psi,  respectively.  The  California  Institute  of  Tech¬ 
nology  Materials  Science  group  reported  a  value  of  185  psi  for  the 
shear  modulus  of  the  particular  sample  tested,  obtained  by  a  torsional 
pendulum  method.  This  is  in  good  agreement  with  the  indentation 
test  results.  Unfortunately,  the  amount  of  material  available  v.'as  not 
sufficient  for  tensile  testing;  therefore,  no  comparison  could  be  rxiade 
as  was  done  with  the  Hysol  specimens. 

A  brief  investigation  was  made  into  the  effect  of  surface 
lubrication  which  was  of  interest  because  the  theoretical  analysis 
assumed  frictionless  contact.  The  results  are  given  in  Table  III, 
Although  the  quantity  of  data  is  rather  meager,  there  definitely  does 
appear  to  be  a  trend  toward  lower  values  of  shear  modulus  with  lubri¬ 
cated  surfaces.  The  data  also  indicate  a  much  greater  difference  with 
the  wet  soap  solution  than  with  the  dry  lubricant  of  talcum  powder 
used  for  the  test  repor  ed  in  Tables  I  and  II,  at  least  for  these  tv.-o 
particular  materials.  The  need  for  s^dditional  investigation  of  sur¬ 
face  conditions  effects  is  clearly  indicated.  However,  surface  lubri¬ 
cation  may  be  inadvisable  with  many  materials  of  interest  because 
of  the  possibility  of  the  material  absorbing  or  reacting  with  the  lubri¬ 
cants  and  thv’.s  undergoing  a  change  in  properties. 


When  the  amount  of  indentation  becomes  an  appreciable 
fraction  of  the  total  specimen  thickness,  the  stresses  and  strains 
produced  may  become  so  large  as  to  cause  specimen  behavior  to 
differ  considerably  from  that  of  a  linear  elastic  solid,  assumed  in 
the  theoretical  analysis.  Tests  were  made  on  several  specimens  of 
different  materials  in  an  attempt  to  determine  the  limit  of  indentation 
for  which  the  analysis  would  still  apply,  as  indicated  by  a  significant 
change  in  the  computed  shear  modulus  from  that  for  lesser  indenta¬ 
tions,  all  other  conditions  being  the  same. 
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table  lii.  ladentatioE  Test  Data  for  Surface  Condition  Study 

For  ail  test  runs: 

isdentor  radius  =0.12 
Penetration  rate  =  2.0 
Teneperaturc  7S''F 

S  in- 
sa.  /min 

Load  readings  taken  after  2  min  relaxation  at  constant  indentation 

Surface 

Indentation 

Load,  P 

Computed 

*  est  Specimen 

Condition 

Depth,  Wg 

(Note  2) 

Shear  Vidilus 

Information 

(Note  1) 

(hi  . ) 

(lb) 

y.  (psi) 

Kysol  8705 

Thickness  =  h  =  0.  504*' 

Dry 

0.097 

12.8 

227 

h/R  =  4.0 

V  =  0. 46 

Soapy 

0.  106 

13.3 

209 

Solithane  1 1 3 

Thickness  =  h  =  0.376" 

Dry 

0.0S8 

11.2 

203 

b/R  =  3.0 

V  =  0.  5 

Soapy 

0.  106 

12.2 

165 

Note  1  -  Dry  =  clean  and  dry,  without  talcum  powder:  soapy  =  coated 
with  wet  soap  solution. 

Note  2  -  Average  of  3  runs  for  each  condition. 

A  definite  limit  applicable  to  all  specimens  was  not  found,  nor 
was  one  expected.  Any  limit  would  depend  on  the  combined  effects 
of  several  variables,  such  as  indentor  radius  and  material  properties 
as  well  as  dimensions,  and  would  likely  be  different  for  every  test 
situation.  The  data  were  nevertheless  helpful,  in  that  they  indicated 
no  significant  deviations  for  penetrations  up  to  one-third  the  specimen 
thickness  for  any  of  the  specimens  tested.  Results  of  deeper  penetra¬ 
tions  were  inconclusive,  but  in  the  case  cf  a  Solithane  113  specimen 
did  tend  to  produce  higher  computed  values  of  shear  modulus.  Deeper 
penetrations  caused  specimens  to  curl  at  the  edges  and  to  assume  a 
cupped  shape,  which  could  be  a  significant  deviation  from  the  geometry 
of  the  analysis.  In  general,  then,  the  maximt.m  allowable  penetra*’on 
depth  will  be  taken  as  one-third  the  specirr.en  thickness,  unless  addi¬ 
tional  data  indicate  otherwise. 
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Following  completion  of  all  the  previously  mentioned  tests, 
the  Hysol  and  Solithane  specimens  were  subjected  to  a  series  of  tests 
to  evaluate  the  effects  of  lateral  dimensions,  particularly  with  re¬ 
gard  to  finding  the  minimum  size  requirea  to  effectively  simulate 
the  infinite  extent  sheet  assumed  in  the  analysis.  Here  again,  no 
exact  limit  can  readily  be  defined  for  the  general  case  because  of 
the  interaction  of  a  number  of  variables. 

The  procedure  consisted  of  repeating  the  same  test  a  number 
of  times,  trimming  the  specimen  to  a  smaller  size  for  each  run.  Data 
from  the  particular  specimens  tested  indicate  no  significant  variations 
in  computed  shear  modulus  so  long  as  the  specimen  length  and  width 
are  at  least  1.  5  times  the  thickness  and  6  times  the  indentor  radius. 
While  this  may  vary  for  other  materials  and  conditions,  there  does 
not  appear  to  be  any  problem  for  the  types  of  specimens  and  con¬ 
ditions  for  which  this  test  device  was  designed. 

2.  Viscoelastic  Materials  Tests 

At  this  time  work  has  only  just  begun  into  the  study  of  visco¬ 
elastic  properties  with  the  indentation  tester  and  very  little  data  are 
available.  However,  one  series  of  tests  has  been  performed  which 
is  of  interest  and  will  be  summarized  in  the  following  paragraphs . 

The  materials  tested  were  different  compositions  of  butyl 
acrylate /acrylic  acid/Unox  221®®  copolymers  formulated  by  Dr. 

A.  R.  Pitochelli  of  these  Laboratories.  Indentation  test  data  are 
given  in  Table  IV.  Note  especially  the  column  of  values  of  Pq/Pj, 

These  figures  compare  the  peak  load,  Pg,  read  immediately  upon 
reaching  the  stated  indentation  depth,  to  the  load,  Pj,  after  two 
minutes  of  relaxation  at  constant  indentation.  The  ratio  Pg/Pz  thus 
is  an  indication  of  ttie  degree  of  relaxation  the  material,  undergoes 
while  under  constant  strain. 

The  same  data  listed  in  Table  IV  were  also  plotted  as  loads 
versus  time  on  log-log  graph  paper,  as  in  Figure  5,  for  example.  The 
results  were  approximately  straight  lines.  Thus  the  load  data  are  quite 
closely  described  by  the  simple  equation 

P(t)  =  Pg  +  at""^  (5) 

^Trademark  of  Union  Carbide  Corp. ,  New  York,  N.  Y. 
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FIGURE  5.  VARIATION  OF  LOAD  WITH  TIME  FOR  INDENTATION  TESTS  OF  BA/AA 
COPOLYMERS 


where  Pg  is  an  equilibrium  load  reached  ^’^ter  a  relatively  long  period 
of  time,  t  is  time,  and  a  and  m  are  constants  determined  for  each  set 
of  test  data  to  give  the  best  curve  fit.  If  the  strain  input  can  be 
approximately  described  by  a  step  function  ^ee  equation  (3)  and  the 
related  discussion.  Section  II. 2  above],  then  a  plot  of  load  as  a  function 
of  time  nrxay  be  directly  converted  to  one  of  shear  relaxation  modulus, 
p.(t),  by  simply  changing  the  scale: 

}i(t)  =  +  bt  ,  (6) 

where  Pg  is  the  equilibrium  mjdulus,  and  b  is  a  new  constamt.  Plots 
of  shear  moduli,  p(t),  versus  time  obtained  in  this  manner  are  shown 
in  Figure  6  for  several  of  the  copolymers  tested. 


Table  IV.  Indentation  Test  Data  for  Various  Compositions 

of  Crosslinked  Butyl  Acr 

ylate /A crylic  Acid  Copolymers 

Composition, 

Shear 

Ratio  of  Peak  Load 

Mole  Ratios  of 

Modulus , 

to  Load  after  2  min 

Remarks  (Applyto 

BA /A A /Unox  221 

psi 

P0/P2 

all  Test  Runs) 

95/5/2. 5 

30 

1.01 

Indentor  sphere 
radius  :  0.  125  in. 

95/5/5 

35 

1.01 

Indentation  rate: 

2. 0  in.  /in.  /min 

90/10/2. 5 

56 

1.01 

Indentation  depth: 
Approx.  0.  10  in. 

90/10/5 

59 

1.02 

Temp:  77'F 

85/15/2.  5 

49 

1.08 

Rel.  humidity: 
Ambient 

85/15/5 

75 

1.  19 

Poisson's  ratio 

85/15/7. 5 

116 

1.20 

of  0.  5  assumed 
for  all  samples 

70/30/2, 5 

40 

4.  18 

Shear  modulus 

70/30/5 

_ 

73 

6.25 

computed  for 
load  after  2  min 

relaxation 

20 


- 


1000 
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FIGURE  6.  SHEAR  RELAXATION  MODULUS  CURVES  FOR  SEVERAL  BA/A  A  COPOLYMERS, 
AS  OBTAINED  FROM  INDENTATION  TEST  DATA 


The  foregoing  indicates  that  the  indentation  test  may  have 
applications  for  investigating  certain  viscoelastic  material  properties, 
in  addition  to  its  present  use  in  measuring  the  equilibrium,  or  elastic, 
shear  moduli.  Properties  which  might  be  studied  include  glassy 
moduli,  relaxation  times,  and  the  slope  of  the  relaxation  curve  in 
the  transition  region.  One  must  remember,  however,  that  several 
assumptions  and  simplifications  were  made,  the  validity  of  which  may 
be  open  to  question.  Any  attempts  to  use  the  method  for  quantitatively 
evaluating  viscoelastic  parameters  must  await  considerable  further 
analysis  and  experimental  confirmation.  Nevertheless,  the  indenta¬ 
tion  test  does  have  possibilities  for  immediate  use  in  qualitative 
studies,  such  as  comparisons  of  the  viscoelastic  properties  of  differ¬ 
ent  materials  or  of  the  effects  of  changing  an  ingredient  in  a  formu¬ 
lation. 


Considering  the  data  of  Table  IV  and  Figure  6  in  the  light  of 
equation  (6)  leads  to  the  following  general  statements: 

(1)  Equilibrium  shear  moduli,  were  found  to  be  directly 
proportional  to  the  number  of  crosslinks  in  the  material,  determined 
by  the  mole  fraction  of  Unox  221.  This  was  an  expected  result, 
based  upon  viscoelastic  theory. 

(2)  The  slope  of  the  relaxation  curves,  p(t),  corresponding 
to  exponent  m  in  equation  (5),  appeared  independent  of  cr  os  slinking 
but  increased  with  increasing  acrylic  acid  content.  The  slope  of 
relaxation  moduli  curves  can  be  expected  to  increase  as  the  test 
temperature  approaches  the  material  glass  transition  temperature. 
Glass  transition  temperature  is  mainly  a  function  of  acrylic  acid  con¬ 
tent  for  these  polymers.  Hence,  it  seems  reasonable  to  conclude  that 
the  acrylic  acid  content,  by  reason  of  its  influence  on  glass  transition 
temperatures,  has  the  greatest  influence  on  relaxation  characteristics 
for  these  copolymers. 

(3)  There  is  indicated  the  possibility  of  superposing  relaxa¬ 
tion  curves  of  samples  of  different  composition  but  having  the  same 
degree  of  crosslinking.  For  example,  the  curve  for  the  85/15/2.5 
composition  could  be  shifted  to  the  right  in  Figure  6  and  merged  with 
the  curve  of  the  70/30/2.  5  composition  to  form  a  single  smooth 
curve.  Thus,  it  would  appear  that  the  shape  and  magnitude  of  the 
relaxation  curve  is  a  function  of  crosslinking  but  its  position  in  time 
a  function  of  acrylic  acid  content.  This  position  in  time  is  indicated 
by  the  value  of  the  coefficient  b  in  equation  (6). 
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Section  V.  CONCLUSIONS 


The  indentation  test  has  been  shown  to  be  a  valid  and  useful 
method  for  determining  the  shear  moduli  of  relatively  low  modulus 
polymers  and  rubbery  materials  in  the  form  of  thin  sheets.  Tests 
conducted  over  a  wide  range  of  conditions  of  indentor  size,  specimen 
size,  penetration  depth,  and  penetration  rate  yielded  data  which,  when 
reduced  according  to  the  theoretical  analysis  presented,  resulted  in 
shear  modulus  values  in  close  agreement  with  those  obtained  by  other 
methods.  The  indentation  test  offers  the  additional  advantages  of 
requiring  only  a  small  quantity  of  material  and  of  the  need  for  little, 
if  any,  specimen  preparation. 

Discrepancies  noted  between  indentation  test  results  and 
moduli  obtained  by  other  methods  can  generally  be  attributed  to 
experimental  variances  and  tolerances.  In  particular,  a  small  error 
in  measurement  of  penetration  depth  can  lead  to  a  much  larger  error 
in  computing  the  shear  modulus,  e.specially  for  cases  involving  a 
large  radius  indentor  at  a  small  depth.  This  probably  accounts  for 
most  of  the  discrepancies.  Another  cause  could  be  friction  between 
the  specimen  surfaces  which  apparently  results  in  slightly  higher  load 
readings  and  computed  moduli.  Also,  the  moduli  from  other  test 
methods,  with  which  the  indentation  test  results  were  compared,  are 
themselves  subject  to  the  experimental  error.  All  discrepancies 
noted  thus  far  are  considered  minor  and  can  be  expected  to  decrease 
as  additional  skill  and  experience  are  gained. 

The  indentation  test  is  a  potentially  useful  tool  in  evaluating 
viscoelastic  properties.  Future  test  efforts  will  be  directed  along 
these  lines,  with  particular  emphasis  on  measurements  of  relaxation 
moduli  and  time-temperature  effects.  Present  plans  call  for.tests 
with  faster  loading  rates  and  longer  durations,  tests  at  various  tem¬ 
peratures,  time- temperature  superposition  studies,  and  continued 
evaluation  of  the  areas  and  limits  of  application  of  the  test  device  and 
theoretical  analysis. 
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A  PPENDIX 

SOLUTION  OF  THE  INDENTATION  PROBLEM 
1 .  The  Lebedev- Ufliand  Solution^ 

Consider  the  state  of  elastic  equilibrium  of  a  linear  elastic 
layer,  infinite  in  extent  but  of  finite  thickness,  resting  upon  a  rigid 
base  and  undergoing  deformation  by  a  rigid  punch.  The  contact  sur¬ 
face  is  a  surface  of  revolution,  and  the  indentation  is  realized  by 
means  of  an  axial  force  (Figure  7).  It  is  further  assumed  that  there 
is  no  friction  between  the  indentor  and  layer,  or  layer  and  base 
(although  the  general  method  will  permit  consideration  of  more  in¬ 
volved  cases,  such  as  those  with  adhesion  between  layer  and  base). 

The  boundary  conditions  imposed  by  these  restrictions  are: 

T  =  0,  w  =  0  when  z  =  h,  (7) 

rz 


T  =0 
rz 


w  =  Wq  -  X{^^)  (j*  ^) )  when  z  =  0, 


a  =0 
z 


(r  >  a) 


wnere  a  and  t  are  the  normal  and  tangential  stress  components, 
respectively,  Wg  is  the  indentation  of  the  center  of  the  punch,  and 
X(r)  is  the  function  describing  the  punch  surface.  Cylindrical 
coordinates  r,z,  Q,  are  employed. 

The  solution  is  facilitated  by  the  use  of  the  Boussinesq- 
Papkovitch-Neuber  displacement  functions.  The  displacements 
are  expressed  in  the  form 

2  fiu  =  -  -  z  ,  (9) 


9$, 


where  u  and  w  are  radial  and  axial  displacements ,  respectively, 
|JL  is  the  shear  modulus,  v  is  Poisson's  ratio,  and  $o  $i  are 
functions  harmonic  in  the  layer  0  <  z  <  h. 

^Lebedev,  N.  N.  and  Ufliand,  la,  S.,  "Axisymmetric  Contact 
Problem  for  an  Elastic  Layer,  "  P.M.M.  ,  U,  3,  320-32b  (1958), 
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FIGURE  7.  GEOMETRY  OF  INDENTATION  PROBLEM 


The  elastic  stresses  are  exp.ressed  in  terms  of  the  functions 
just  introduced  as 


a  =  2(1 
z 

-'■>  'sr 

8^», 

■ 

(11) 

and 

d 

"^rz  dr 

-‘?rj  • 

(12) 

Substituting  (7)  and  (8)  into  the  relations  (9)  through  (12)  gives 
the  boundary  conditions  in  the  following  form: 


The  harmonic  functions  ^>o  ®^re  taken  in  the  form 


where  A(X)  and  B(X)  are  functions  to  be  determined  and  Jo(Xr)  is  a 
Bessel  function. 
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^0  <  t  < 


(23) 


f  jyXr)  sis  >x  d> 
'^0 
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, 
t  S 


r*)  ^  tc  >  r) 


cos  M  d>. 


(0  <  E  <  r) 

f:  >  r) 


? 
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COS  (fXrsin  c5 


(26) 


(27) 


Upoa  istegraEisg  tie  righs-hatnc  portion  of  JZ3)  by  parts.  substitnEing 
the  resalting  expression  for  A(X)  utto  ecoation  (22),  and  making  ase  of 
formsila  (23).  •we  may  sbov  that  equation  (22)  is  satisfied  identicaliy- 
Detaiis  of  these  steps  wiii  not  be  shown  here. 

Sebstitudon  of  Af>. )  into  equation  (2!)  and  making  use  of  the 
formulas  {!c)  and  (27)  leads  to  the  reiadon 


g(X)  cos  Xt  cos  {Xrsin0)  dX 


1-v 


ho  -  X(r)J 


(28) 


Introducing  a  new  integration  variable  t  =  rsin  9  into  the  first  integral 
and  ntanipulating  (28)  further,  there  results 


(Ji(rsin  9)  -  — 


{ G(t  rsin  9) 


G{t-rsin  0)] 


Hr) 


(29) 
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wnere 


C-Cx)  =  f 


gfX  )  cos  Xx  d>-  . 


tise  Fc^irier  cosine  transform  of  the  function  g(X),  Putting 


-  T"  ® 


F(x)  =  6«x)  -  —  \  6(e)  (G(t  -  x)  -  G(t  -  x))dt  .  (32) 

then  equation  (29)  wili  have  the  form  of  Schlomilch’ s  integral  equation^ 


^  F(rsinf)d6  =  f(r) 


which  has  the  solution 


F(x)  =  ^  f 


c 

(0)  *  X  \  f'(xsm0)d6  , 
*^0 


where  the  prime  denotes  the  derivative.  The  result  of  the  foregoing 
is  an  integral  equation  for  the  unkno;»-n  function  6(t): 


-  ~  \  <i>(t)  [  G(t  -  x)  -r  G(t  -  x)  ]  dt 

'^0 


2^  r  c'‘  1 

=  j  ^  j  wg  -  X  J  X'  (x  sin0)  d^  ,  0  <  x  <  a  .  (35) 


Vhittaker,  E.  T.  and  V/atson,  G.  N.  ,  A  COURSE  OF  MODERN 
ANALYSIS,  Fourth  Edition,  Cambridge  Un  versity  Press,  229  (1958). 
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if  it  is  possible  to  find  a  solution  to  the  integral  equation  (33)  in  the 
form  of  a  function  with  a  continuous  derivative,  then  equations  (18). 
(19),  (20),  and  (23)  give  a  complete  solution  of  the  contact  problem 
under  ccH^sideration. 


iSasy  quantities  of  interest  can  be  expressed  immediately 
in  terms  of  the  function  ^it).  omitting  the  intermediate  forrrmlas. 

For  example,  taking  the  general  expression  for  stress  as  given  in 
(11),  using  the  functions  defined  in  (13)  and  (19).  and  making  the 
substitutions  of  (20),  123),  and  (24),  gives  the  expression  for  the  dis¬ 
tribution  of  normal  stresses  directly  beneath  the  punch,  on  the  sur¬ 
face  z  =  0: 


(a  ] 

'  z‘z=0 


^  <y(t>dt 


d>(a; 

Va-  -  r" 


r  <  a 


(36) 


integrating  (36j  over  the  area  of  the  circle  of  radius  a,  we  find 


C 


P  =  2t  \  d>(t)  dt 


(37) 


in  which  P  is  the  total  applied  force  to  the  punch. 
2-  AoDlication  to  the  Indentation  Test 


In  an  indentation  test  the  only  easily  measured  quantities 
are  the  indentor  dimensions,  the  maximum  penetration  depth  w^,, 
and  the  applied  load  P.  Therefore,  to  be  of  any  use  to  the  experi¬ 
menter  the  Lebedev-Ufliand  solution  must  be  in  a  form  which  permits 
direct  calculation  of  the  shear  modulus  from  these  measurable 
quantities . 

The  following  analysis  is  based  upon  the  use  of  a  spherical 
indentor,  which  has  the  contour 

X(r)  =  R  -  Vr'  -  ?  .  (38) 

where  R  is  the  radius  of  the  sphere.  Then  equation  (35)  takes  the 
form 
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From  the  condition  of  continuity  of  stresses  at  radius  a,  it  can  be 
deduced  that 

d(a)  =  0  =  ^l)  .  (43) 

Applyiijg  this  to  equation  (40)  yields  the  expression 


Substitution  of  (44)  into  (40)  gives 


which  is  a  Fredholm  integral  equation  of  the  second  kind. 

The  problem  now  becomes  one  of  solving  the  integral  equation 
(45)  for  the  function  tj(|).  Once  this  has  been  done,  the  dimensionless 
load,  (l-v)P/uR^,  and  dimensionless  deflection,  Wq/R,  can  be  ob¬ 
tained.  In  the  general  case  it  is  necessary  to  use  numerical  methods. 


Numerical  Solution  of  Equations 


Solution  of  equation  (45)  involves  first  finding  values  for  the 
function  K(u)  in  the  interval  0  <  u  <  2.  The  function  K(u),  equation 
(41),  itself  involves  the  integral  of  an  oscillating  function  over  an 
infinite  interval  which  cannot  be  evaluated  in  closed  form.  Filon's 
method  was  found  to  perform  the  indicated  quadrature  adequately  by 
evaluating  K(u)  for  the  finite  interval  0  S  a  s  10  in  increments  of 
0.2.  This  was  programmed  as  a  subroutine.  For  the  solution  of 
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the  integral  equation  (45),  a  Gauss- Legendre  quadrature  was  used. 
In  this  manner  (45)  may  be  rewritten 


A  Mr.)  [K(t.  -r 
J  J  J 


l)r  K{t.  -1)1 


(46) 

r  I 

In 

!-r  pT. 

1 

1-pT. 

1/ 

in 

where  the  A  -  are  weighting  coefficients  and  the  and  tj  are  roots  of 
the  Legendre  polynomial  of  degree  N.  The  resulting  set  ot  equations 
is  solved  for  values  of  u  at  |  =  These  values  of  u  are,  in  turn, 
used  in  evaluating  (42)  and  (44)  by  Gauss -Legendre  quadrature.  The 
results  are  presented  in  the  body  of  the  report  (see  Figures  2  and  3) 
as  curves  of  dimensionless  load  (1-v)P/hII^  versus  dimensionless 
indentation  Wq/R  for  various  values  of  h/R,  the  ratio  of  indentor 
radius  to  slab  thickness. 
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